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ABSTRACT 


This report concludes a multiphase program to demonstrate the utility of 
remote sensing data analysis to water quality monitoring. This phase demonstrates 
that data analysis and image processing techniques can be applied to Landsat 
remote sensing data to produce an effective operational tool for lake water 
quality surveying and monitoring, The major findings of this report area 

(1) Digital Image processing and analysis techniques have been designed, 
developed, tested and applied to Landsat Multispectral Scanner (MSS) 
data and conventional surface acquired data. Utilization of these 
techniques can facilitate the surveying and monitoring of large numbers 
of lakes in an operational manner. 

(2) Supervised multispectral classification, when used in conjunction with 
surface acquired water quality Indicators, can he used to characterize 
water body trophic status. 

(3) Unsupervised multispectral, classification, when interpreted by lake 
scientists taaillar with a specific water body, can yield classifica- 
tions of equal validity with supervised methods and in a more cost- 
effective manner. 

(4) Image data base technology can be used to great advantage in character- 
izing other contributing effects to water quality. These effects 
include drainage basin configuration, terrain slope, soil, precipitation 
and land cover characteristics. 


INTRODUCTION 


This publication present b the results of the final work completed through the 
joint efforts of JPI.'s Image Processing Laboratory (I PL) and U.S. Environmental 
Projection Agency, Las Vegas (EPALV) on the Lake Glassification Task. The report 
lias been divided into two major sections. 

Section I, entitled Significant Technology, presents brief descriptions of 
four digital image processing software programs on systems (LAKELOC, STORE!', 

STATS 2, and IBIS). LAKELOC, STORET and STATS 2 were" developed specifically for 
use by the lake classification projects. The IBIS system was developed originally 
for land use planning activities but subsequently successfully applied to the lake 
classification effort. These four systems formed the technical core through 
which all Landsat data and much control sensed data were analyzed. 


Section II, entitled Technology Applications, presents detailed discussions 
of three major lake classification studies which were conducted between 1977 and 
1979. The studies are presented separately in chronological order beginning with 
the Lake Mead Intensive. Area Study. This study was designed as a demonstration of 
the application of Landsat data to water quality monitoring of a single, large 
water body. Multlspectral classification was performed and the results compared 
with conventionally acquired contact-sensed measurements gathered on corresponding 
Landsat flyover date/- 1 . The Lake Tahoe Study was conducted as a demonstration of 
the application of IBIS to a water quality monitoring situation. A data base was 
constructed consisting of data elements from Landsat, digital terrain data and 
Tahoe Regional Planning Agency statistics. Data base manipulation and some pre- 
liminary modelling was performed for demonstration purposes. The third study, 
involving lakes In Montana, Utah and Michigan, was the final phase of the Lake 
Classification Project designed as a demonstration of computer program capabilities. 
Multlspectral classification was performed on some lakes selected for this study. 
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SIGNIFICANT TECHNOLOGY 


Tlnee ho! t wart* systems were designed specifically to handle* the tank of lake 
classification, LAKELOG and its associated follow-on programs FARINA and STATUS 
were developed to aid the analyst with Identification and Isolation of lakes 
selected for study. Either Landrat or digital aircraft Imagery could be used as 
input to LAKELOG , The two follow-on programs could be used to extract information 
for each lake of interest from each spectral band and present tills information in 
the form of a computer listing. 

The STORKT system was designed to present contact-sensed water quality mea- 
surements in the form of scatter diagrams. The Information source for this 
program was data from the National Eutrophication .Survey, 

The computer program STATS 2 was developed to be used interactively to locate 
and identify sample, site locations on lakes selected for study. Multispectral In- 
formation pertaining to each sample site could be printed out. In addition, super- 
vised classification could also be performed using 8TATS2 by using sample site 
statistics as input to an automatic classifier. 

The IBIS system was developed at ,JPL originally for land use studies requiring 
the development of information data bases to be constructed from varied sources. 
IBIS was used in the lake classification project to build a data base composed of 
information from sources such as land sat, Defense Mapping Agency digital terrain 
data and maps provided by the Tahoe Regional Planning Agency, 


A. LAKELOG 

As described in previous project reports (Refs. 1-1, 2., 3, b ) , lake Isolat ion 
performed at the XPL relied solely upon batch processing techniques. Analysts 
used only the IR2 band from Landsat imagery to define water by examining histograms 
and setting an average threshold (Ref. 1-5). This method, while fairly accurate, 
for defining water boundaries and aiding the lake extraction process, was inherent- 
ly slow. The analyst was required to wait for processing of hardcopy photo 
products before continuing with each phase of lake extraction. In addition, small 
lakes were extremely difficult to separate, especially if these lakes were directly 
associated with river systems or closely situated to other water bodies. 


To streamline and increase the accuracy of lake extraction, analysts associ- 
ated; with this project developed an interactive system of lake Isolation. The 
system consisted of three interactive programs which reduced the time needed to 
locate and isolate a lake, increased the accuracy and sensitivity of the water 
detection scheme and output a statistical and surface area listing for each lake 
to be surveyed. The system was called LAKELOG. The following section describes 
the hardware configuration, program operation, water detection algorithm and 
hardcopy output associated with the LAKELOG system. 
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Hardware 


The host computer is currently an IBM 370/158. The display controller 
used is a Ramtek G100B, a versatile video display device that can be used to dis- 
play gray level images and graphics data. The Ramtek is a solid state refresh 
memory system with a display format of 512 lines by 640 elements. Readback from 
the refresh memory is available under software control. It is possible to display 
6-bit gray level Images along with two graphics planes, and the user may selective- 
ly write or erase the displays point by point. Manipulation o i. the graphics data 
can be accomplished with the aid of a trackball cursor. Figure 1-1 illustrates 
the configuration of the Interactive hardware as it is arranged for the operati a 
of the LAKELOC program. 


2. Operation of LAKELOC 

For the purposes of lllustrat on, we will follow the operation of LAKELOC 
as it would be applied to a scene in southern Illinois. Although any number of 
lakes could be extracted fr m a scene, we will limit the operation to one lake, 
in this case, Crab Orchard Lake. 

For a given digital data scene, such as Landsat, the user may selectively 
display 5.12 y 640 element subareas until he locates the water body of interest. 
Automatic linear contrast stretching of the displayed scene can be performed 
during this operation to aid in the location of the lake. 



Figure 1-1. User console with video display device and track ball 
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Once a lake has been located, the trackball cursor Is then set on the desired 
lake and a default SO x SO element box Is drawn on the graphics plane about the 
cursor position. Figure 1-2 Illustrates the default box drawn about the cursor 
pocltlon on Crab Orchard Lake. Since only the area within the box will be acted 
on by the watci detector, the user must correct the size and position of the box 
relative to the lake so that the lake Is contained within the box boundaries. The 
size Is ctianged by a simple command to the program, which allows the manipulation 
of the trackball cursor to control the box dimensions. The position of the box 
Is also controlled In the same manner by the trackball. An automatic mode can 

also be used In which the user places the cursor on the lake of Interest and gives 

the command "window". A box Is automatically drawn about the lake boundaries 
with no need for trackball manipulation. If a lake encompasses a major portion of 
the video display, a command can also be given to the program to "window" the en- 
tire screen In an automatic mode. This allows the user to forego manipulation of 
the cursor if the entire screen dimension is desired. Figute 1-3 illustrates 
Crab Orchard Lake completely enclosed by the box after manipulation by the cursor. 

Once the box has been satisfactor y positioned about the lake boundaries the 

water detector can be used to Isolate tue water body in a binary form. (a 

detailed description of the water detector follows In the section entitled "Water 
Detection Algorithm". ) In the binary form, the water bodies appear as white and 
non-water features as black. The user can also magnify the area within the box 
boundaries by Issuing a "zoom" command with the appropriate magnification factor. 



Figure 1-2. Crab Orchard Lake with default 30 x 50 
element box 
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Figure 1-3. Crab Orchard Lake with correctly 
sized and positioned box 


rhe zoom command thus redisplays a magnified picture of the boxed area only, 
directly over the existing image. At tnj conclusion of the editing session the 
nuignit ied image »s erased from the screen, leaving the original image. Ibis 
allows the user to continue uninterrupted with other lakes contained in the ex- 
isting scene. 


Magnification of the image allows the user to easily determine the exact 
boundaries of the lake as opposed to any extraneous water information which may 
also be displayed in the scene. The task of editing out extraneous pixels has 
been in t he past the most time-consuming cho e in the water analysis project. 

Aided only infrequently by a map, the user had to decide what constituted the 
boundaries ot the lake in question. Previously the user relied on pixel listings 
and hardcopy photographs to locate the lake boundaries. In the case of a large 
lake, one was often hampered by cumbersome pixel listings which had to be care- 
fully pieced together to recreate a lake image. With LAKE LOG , the magnification 
factor in conjunction with the easily manipulated trackball allowed the user to 
perform the editing task in a matter of minutes as opposed to a duration of 
several days. Figure 1-4 illustrates Crab Orchard Lake in binary format. The 
detached white areas represent extraneous water information which is not associated 
with the lake. 
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Figure 1-4. Crab Orchard !«»ke In binary form with magnification factor 
of 2 , before editing of extraneous water features 


The removal of water bodies not associated with the lake ol interest can be 
performed in three ways. In the first method the t rackball-cont rolled cursor is 
set point -by-point on the areas to be removed. The default size of the area re- 
moved is one pixel; however, the user can specify the number of surrounding pixels 
to be removed for each erase operation. This method is most useful when working 
in close proximity to the boundaries ot the lake of interest, when it is most 
imperative not to remove too large a section of pixels. The second method is 
used for removing closely spaced extraneous pixels which are not in close prox- 
imity to the lake. This method utilizes continuous erasures as the trackball 
cursor is moved across the screen. The size of the area cbout the cursor position 
to be erased can also be controlled by the user in this mode. The third method 
allows the user to place the cursor on the lake of interest and isolate it auto- 
matically. The user must first, if necessary, "detach" the lake from any other 
water bodies by one of the two previous methods. A command is then entered which 
quickly erases all water information contained within the boxed area except for 
the lake on which the cursor is situated. Figure 1-5 depicts Crab Orchard Lake 
still in binary format after all extraneous information has been removed during 
the editing phase. 

Once the user is satisfied that he has isolated the lake of interest, he then 
assigns the lake a name and commands the program to save the binary image of the 
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Figure 1-5. Crab Orchard l.ake in final edited form 


lake on a disk data set. The Lake's position in the disk data set is exactly 
the same as it is in the original Landsat scene. LAKELOC returns to the user the 
exact position and size of the extracted lake, image as it appears on the disk 
data set, and a parameter data set is also created which contains this positional 
information. This information is necessary to the operation of the follow-on 
programs for LAKELOC. At this time the user is able to continue processing any 
number of lakes or, if finished, he can fetch the follow-on programs which will 
process the statistical data. 


3. Follow-on Programs 

The output from LAKELOC consists of a binary mask disk data set con- 
taining the extracted lakes and a parameter data set containing the positional 
information and lake names. The output size of the binary disk data set is 
exactly the same as the size of the original landsat image used as input to 
Lakeloc. in the next step, the binary data set is used by the program FARINA to 
make out of each corresponding spectral channel in the original I^andsat frame 
each water feature which has been processed through LAKELOC. The output is four 
data sets containing the original l)N values for each lake in each of the spectral 
bands. This output can in turn be used as input to the program STATUS and as in- 
put to follow-on MSS classification programs. 
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STATUS, utilizing the parameter data set from LAKELOC or punched parameter 
cards, produces a statistical analysis in the form of a hardcopy printer listing of 
lake statistics in all four spectral channels. The lakes are listed by name and 
ranked according to size. Two tables are printed, the first consisting of lake 
statistics such as pixel count, surface area calculations and shoreline perimeter 
calculations. The second table lists lake MSS statistics, such as the moan DN 
level for each lake in each spectral channel and the corresponding standard 
deviations. Figure 1-b is a reproduction of the lake statistics table, including 
the statistics for Crab Orchard lake among other lakes processed from the same 
Landsat frame. Figure 1-7 reproduces the lake MSS statistics for the same lakes 
listed in Figure 1-b. 
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4. Water Detection Algor itlun 

In the paHt, the detection of pixels whose Instantaneous field of view (1F0V) 
Is ttiat of water wus accomplished by threshold lng handset's Infrared band 7 (0.8 - 
1.1). The low reflectance of water In this spectral range conveniently produced 
a blmodal distribution of DN's: one peak for water, another peak for non-water. 

This technique was accurate except In the case where the 1F0V of the scanner was 
viewing a combination of water and non-water areas such as at the shoreline of a 
lake. In this situation, the problem became one of trying to estimate the pro- 
portion of each material In the IFOV. 

Horowitz (Ref. 1-b) Work and Gilmer et al. (Ref. 1-7) have investigated the 
proportion estimation problem and have encouraging results. Work and Cllmer have 
estimated the proportions of water, bare soils and green vegetation using Landsat 
bands 5 and 7. This technique requires an estimate of the spectral signature for 
pure water, pure bare soil, and pure vegetation. While the spectral signature of 
water is fairly easy to estimate, that for soil and vegetation becomes more 
difficult. The many variables involved, such as different soil types, vegetation 
cover types and thickness of the vegetative cover cause considerable error when 
estimation is attempted by a completely automatic processor. 

An alternate approach, and the one chosen to implement in the LAKELOC system, 
considers the mixture classes to be only water and non-water. Bands 5 and 7 are 
used In the detection process, as it was found that bands 4 and 6 offer little in 
additional information. The estimation of the spectral signature for water and 
non-water Is made over a region within and immediately surrounding the water body. 

The spectral signatures (mean DN values) were estimated by an iterative pro- 
cedure. First the 2-d imenslonal space (band 5 vs. band 7) is partitioned into two 
regions in which the populations of water and non- *ater typically cluster and the 
mean is then recomputed for those DN's which fall within a neighborhood of the 
Initial mean. This process is continued until a convergent mean has been found for 
each region. 

The proportion estimation Implemented uses a technique proposed by McCloy 
(Ref. 1-8). In Figure 1-8, W is the mean for water, U is the mean for non-water 
and P is the DN for any given pixel. P' is the projection of P onto the line 
segment WU. If /WU/ is the length of the line segment WU and /WP '/ is the length 
on line segment WP", then the proportion estimate q for water is: 


<1 


/WP V 
/WU/ 


where 0 * q 1 . If P' does 
of the closest point W or U. 
pixel is defined to be water 


not fall between W and U 
A decision threshold is 
or non -water. 


it is given the position 
set for q at which the 


B. STORET 

The STORET work was begun by JPL in late 1978 and left uncompleted due to 
termination oi the project. The work was attempting to develop experimental 
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BAND 5 

Figure 1-8. Geometric Interpretation of water 
detection algorltlun 

models that could relate remotely sensed mu 1 1 i spec t ra l data to actual water 
quality measurements. 

KPA produced for JPI. a set of STORET data that contained the chemical com- 
position of water samples along with Secchi depths and turbidity measures for 
many of the fresh water lakes in the U.S. for the years 1973, 1974 and 1975. 

Most sample sites had samples taken during summer, fail and spring months. 

Computer software was developed that would allow the analyst to extract from the 
STORET files specific water quality measurements for a given year and season. 

Due to the high volume of data, an analyst would typically pick some subset of 
the 21 possible water quality measures for analysis. 

In order to find relationships between remotely sensed data and water quality 
measurements, some basic software tools had to be developed. A data format 
structure was defined that would allow both water quality measures and remotely 
sensed data to be stored concurrently in one data file. Each measure was con- 
sidered as a separate variate. The name of each variate was stored in a header 
record for the file. This scheme allows much flexibility in both the data format 
and the software that manipulates the data. 

1. basic Software 

Some basic software developed for analysis of this multivariate data is 
now described. 

(1) Plots. A program that could plot, on CRT graphics plane, scattergrams 
of selected pairs of variates. For example, chlorophyl-A might be 
plotted against total organic nitrogen; Secchi depth vs l«indsat MSS-5, 
etc. These plots could be made in linear or log domain, with manual 
or automatic scaling of axes. A linear regression line could also be 
plotted along with the correlation coefficient of the plotted points. 
Provisions existed for easily producing hard copy prints of scatter- 
grams after they were plotted on the CRT (see Figures 1-9, 10, 11, 12). 
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Figure 1-9. A scattergram of the log of total organic nitrogen (TON) 
on the X-axis versus the log of chlorophyll a (CHLFA) on 
the Y-axis. Data is from STORF.T for the summer of 1973. 
The correlation coefficient is also given 
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Figure 1-10. A scattergram of the log of total organic nitrogen (TON) 
versus chlorophyll a (CHLFA) . Included with the plot 
are the linear regression line and the slope. Data is 
for the summers of 1973, 1974, 1975 
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Figure 1-11. A scattergram of the log of total organic nitrogen (TON) 
versus total phosphorus (PTOT) for the summers of 197^, 

1 974, 1 975 
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I iRurt* 1-12. A scattergram of the log of total phosphorus (PTOT) versus 
ortho-phosphorus (POKTHO). Note the expectedly high 
correlation coefficient 
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(2) Thresliold ing . It was often desirable to disregard samples which had a 
certain variate above or below some threshold value. Since turbid water 
lias greatly different reflectance character 1st ics tlian nonturbld water, 
it was common practice to discard those samples with turbidity above 
some level. 

(3) Normalization. Most water quality measures are In different units, e.g., 
mg/1, meters, etc.; correlation between them Is difficult to ascertain. 
Normalization of each variate based upon Its mean value and standard 
deviation allowed relationships between variates to be observed more 
read lly . 

(4) Ratio ing. In many cases it Is the ratio between water quality measures or 
reflectances that is Important. A program was developed that allowed new 
varlater to be computed that wore arithmetic functions of existing variates. 
An example of a useful ratio would be the ratio between total organic ni- 
trogen and total phosphorus. 

(5) Principal Component Transformation. Correlation between many variates 
could be analyzed by performing a principal component transformation on 
selected sets of variates. This technique was also being investigated as 
a possible method of removing noise in the multivariate data. 

2. Work Plan 

At the time that the STORET work stopped, JPL. liad progressed to a point of 
trying to understand the interrelationships of the water quality measurements 
themselves. The aforementioned software had Just been completed and an analyst 
was being trained in its use. The work plan for the future was as follows: 

(1) Develop a reliable trophic state index from the water quality measure- 
ments. This implies the ability to compute a number whose value 
describes the trophic state of the water body from which the sample 
was taken. The STORET data gave us an excellent data base to work from 
and it is imagined would allow development of a trophic state index that 
could be applied over a larger regional area. 

(2) Predict the trophic state index from remotely sensed data. Cnee the 
trophic state index was computed, it could be placed into our data base 
as an additional variate; then the existing software could be used to 
correlate mult ispectral remotely sensed data to the trophic state index. 

(3) If the trophic state index could be reliably predicted, the next step 
was to predict the water quality measures from the remotely sensed data. 
More realistically, ratios or other relationships of the water quality 
measures are probably the best that could be done. 


C. STATS2 

During the course of the water quality task, it was frequently necessary to 
compute training area statistics of the mult ispectral imagery at precise locations 
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within the scene. These lorat (unit were typically those at which EPA personnel 
had taken water samples during the satellite or aircraft pass. 

To facilitate the training area vocation procaaa, an Interactive statistical 
training area program wan developed called STATS2. With STATS2, the analyst was 
able to view the scene on a CRT, Interactively choose the spectral channel that 
gave the beat rendition and cont rant -enliancr the scene. Through the use of s 
t rac kba 1 1 -cont rol led cursor, the analyst could then outline the precise location 
tor which mult lspuctral statist lea were desired. The program would then compute 
the statistics and save them on a dlsx file. 

3TATS2 also contained capabilities for editing, updating and merging of 
statistics. in addition, the analyst could perform a "first-cut" mult lspectral 
c lasslf Icat Ion ao as to decide whether or not the training area was well chosen 
and would perform well in a later r lass 1 f lcai f on procedure. 


I). IBIS 

Although the Image Baaed Information System (IBIS) was not designed specif- 
ically for the water quality monitoring tasks, one of Its major Initial applica- 
tions was directly connected with watershed monitoring in the l^ke Tahoe Basin. 
The IBIS system was originally designed as u computer-based approach to the 
anal><«ih of geographical situations such as those associated with land use plan- 
n ing act ivlt lea . 

L. System Description 

IBIS la composed of general purpose and specialized computer programs which 
can be grouped into logical steps to build an information data base. Functional- 
ly, IBIS represents a selection of programs which operute under the IPL'a VICAR 
Image Processing System. The IBIS system la raster (image) based. Most data 
entered into thv. system are ir. •"aster format. However, the system Is also de- 
sigied to allow Integration of graphical and tabular data aa well (Ref. 1-13). 

logical and mathematical interfaces arc provided in the system to link all 
data files In an IBIS data ba-r superstructure. Figure 1-13 Is a configuration 
diagram depicting the IBIS data base concept. By manipula' ing the Interfaces, 
information can be derived from simple associations to or comparisons between two 
or more data flies stored In the data base. More complex procedures, such an 
polygon overlay and cross-tabulations can be performed to acquire other signifi- 
cant data. 

The Image formatted data plane Is the primary data type utilized in IBIS 
processing. IBIS data planes are obtained directly In Image form, ns from Land- 
sat Imagery, or the planes can be derived from data compiled by sources such as 
the U.S. Department of Agriculture, U.S. Bureau of the Census and the Defense 
Mapping Agency. 

Most image data sets entered into the data base are derived from Iwindsat or 
other mult lspectral scanner sources. Since image data planes can be derived from 
different sources, often no common spatial alignment exists between them. The 
system thus provides a means of registering these data sets to obtain a unified 
spat ial surt ace. 
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Figure 1-13. A configuration diagram (from Ref. l-9) 
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Graphical or vector data can aiao be entered Into the IBIS data bane. Graph- 
ical data, hucIi as those obtained from conventional mapa, are electronically 
dixit ized on a coordinate digit izer. Ah with image data, the graphical data files 
must be in registry with the primary data base. Geometric corrections are per- 
formed before the graphical files are transformed Into image space uclng a two- 
step process. First, a general surface fit is achieved through the use of a 
least squares affine transformation. Then an exact geometric correspondence to 
the primary data base is obtained through a "rubber sheeting" procedure. Defor- 
mation of the original surface is controlled by selection of tlepolnts which link 
geographical features identifiable on both the graphical data file and the 
primary data base. 

Tabular data can be entered into IBIS through computer cards or digital tape. 
These data are stored in a tabular file which is linked to the data base through 
a logical Interface. 

The most Important element of the IBIS concept is the gem cference plane. 

The georeference plane is a map-based graphical representation of areas ot 
interest, such as drainage basins in a watershed area. The georeference plane 
is constructed In the form of a polygon file which is registered to the primary 
data base and used in data aggregation and map generation procedures. 

Once the georeference plane is transformed into image space, each polygon or 
region of Interest is identified by assigning a unique data number (DN) tc each 
region. Alter region identification, the geeret erence plane is used in higher 
order IBIS procedures such as polygon overlay or modeling. Each Image plane 
is referenced to one or more georeference planes. All tabular files are also 
linked to at least one georeference plane. By its unique construction, the IBIS 
data base provides the user the ability to manipulate data from several sources, 
which despite their original disparity, are all referenced to a common base. 

Figure 1-14 provides a conceptualization of the registered data plane format of 
the IBIS system. 


2 . Application of IBIS to Water Quality Studies 

Major emphasis in past water quality studies conducted at the I PL had been 
placed on assessing the viability of l.andsat as a monitoring device. A useable 
system for lake extraction and classification was developed as a result of these 
tasks. However, lakes were only examined through the comparison of certain water 
quality measurements with mull lspectral scanner data analysis. As a result, not 
much emphasis was placed on locating and defining any contributing factors, such 
as the Impact of non-point source pollutants ftom the surrounding land mass. 
Possible sources for pollution were suggested only for the purposes of discussion 
in the analysis of classification results. 

In order to conduct a more comprehensive study which attempted to define the 
causes and dynamics of lake pollution using remotely sensed data, the resources 
of IBIS were essential. Remotely sensed data, such as l*andsat imagery, required 
augmentation by conventionally acquired data, such as hand-drawn maps, runoff 
measurements and precipitation data, if an accurate data base were to be produced. 
Similarly, conventionally acquired data were greatly entianced when integrated with 
spectral Information provided by satellite. The task of combining these data, 
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Figure 1-14. Conceptual diagram depicting data planes in registration forming 
the IBIS data base 
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however, was enormous unless a system were available which could tie all data 
elements to a common base. Rather than viewing a precipitation map aide by aide 
with a Landaat scene, it was more desirable to examine the map as it appeared 
"overlayed" on the Landsat image. This technique could be easily accomplished 
using the IBIS resources. 

Tremendous amounts of environmental data have been gathered for the water- 
sheds of our nation's larger lakes. In the case of the Lake Tahoe Basin, the U.S. 
Department of Agriculture had conducted an intensive study of the basin's water- 
shed in an effort to define controls on future land use. These data provided an 
excellent source of material which could be integrated with Landsat MSS and digital 
topographic data to not only test the IEIS concept but to add a new dimension to 
the water quality management projects. 
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SECTION II 


TECHNOLOGY APPLICATIONS 

This section details the techniques used and significant results of three 
demonstration projects conducted during the period of 1977 through 1979. The 
lake Mead intensive Area Study was completed In 1978. This study applied 
digital Image processing techniques to Landsat data acquired over the Lake Mead 
area. Mult ispectral classification was performed and the results were compared 
with contact-sensed water quality measurements gathered by EPA personnel. 

The lake Tahoe study applied digital Image processing techniques of the 
IBIS system to build an information data base for the Lake Tahoe watershed area. 
The data base was used to demonstrate the capabilities of the IBIS system and Its 
applicability to a watershed monitoring situation. I^ndsat based Images were 
created which combined conventional maps with Landsat data. Additionally some 
preliminary modelling of terrain data was attempted. 

The study of selected lakes in Montana, Utah and Michigan was undertaken as 
a final demonstration of digital lake isolation and classification techniques 
using Landsat data. Mult ispectra 1 classification was performed on some lakes 
selected for this study. Results are presented only in terms of digital image 
processing products as no analysis has yet been made by EPA personnel,. 


A. LAKE MEAD INTENSIVE AREA STUDY 
1. Purpose of the lake Mead Study 

The study of lake Mead was performed in order to establish the feasibility of 
monitoring the trophic condition of a water body over a specific period of time. 

This project was based upon the application of landsat MSS data and airborne sensor- 
acquired data. These remotely sensed data were supported by ground truth measure- 
ments gathered by EPA l iinno legist s . The period of landsat coverage for the study 
was from June through September L97b on lake Mead, Nevada. 

Trophic classification of lake Mead proved to be most successful when performed 
by an unsupervised clustering routine. The results of unsupervised classification 
and a comparison with ground truth measurements are presented within this section. 

The promising results obtained from the Colorado Lakes Study (Ref. 1-2) and the 
Illinois lakes Task (Ref. 1-3) led to the definition of a project in which a single 
water body of considerable size would be more closely monitored. The intent of the 
project was to apply the skills and techniques which had been developed by the IPL 
during prior lake classification investigations to a new and unique situation. In 
the past the lakes surveyed had been extremely small, sometimes limited to only 30 
pixels in a landsat frame. Much effort had been expended in the development of the 
image processing skills necessary to perform the task of lake identification, extrac- 
tion and analysis in a timely manner (Ref. 1-5). The 'ake Mead Intensive Area Study 
offered a test of the techniques which liad evolved from the previous work. In ad- 
dition, the study was oriented as a temporal investigation in which one relatively 
large lake would be cliarac ter ized through the correlation of ground truth measure- 
ments and landsat MSS analysis, and classified as to both trophic classes and tur- 
bidity Levels on a space-time basis. 
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2. The Lake Mead Study Area 

lwike Mead lb the largest man-made lake In the United States. It was created 
by the Interstate Impoundment of the Colorado River by Hoover Dam In 1935. Lake 
Mead, the reservoir formed by this dam, was completely filled in by 1941 (Ref. 
11-1,2). Located in the Mojave Desert, as seen in Fig. 11-1, the lake bed con- 
sists mainly of sand, silt and clay Joined with some expansive areas of soluble 
beds of gypsum and reck salt (Ref. 11-2). The multiple use of water from Lake 
Mead results in widely fluctuating water levels, which generally follow a pattern 
of slow drawdown from midsummer to the following spring when a rapid rise occurs 
(Refs. 11-3,4). The water is hard, with total dissolved solids of about 800 ppm. 
The annual temperature cycle of Lake Mead can be classified as warm monomlctlc. 

This is characteristic of lakes of he warmer latitudes in which the temperature 
of the water never falls below 4°C at any depth, one circulation occurs each year 
in winter and the lake is directly str .if led in summer (Ref. 11-5). The 
reservoir provides municipal and Industrial water, water for irrigation, hydro- 
electric power and a recreational area visited by over one million people each 
year (Ref . 1 1-6) . 

Due to the proximity of Lake Mead to the USEI’A's Environmental Monitoring 
and Support Laboratory at las Vegas (EPA/LV), sampling surveys could be convenient- 
ly and cost-effectively conducted. Further, lake Mead was relatively large when 
compared with lakes previously studied and meteorological conditions of the 
lake Mead area seemed favorable in terms of cloud cover and ice conditions, which 
can severely limit the extent of Iandsat scene availability. One of the most 
Important factors determining the selection of Lake Mead was that the reservoir 
is situated within the 14% image overlap of the Iandsat sensor scan path. This 
allows coverage on successive days, thereby providing not only more images to 
select from, but also an opportunity to test the validity of trophic classifica- 
tion from one day to the next. 


3. Sampling Procedures 

following the selection of lake Mead as the study site, the USEPA provided 
sampling crews to measure water quality parameters on the lake from early June 
through September 1976. Samples were collected by Bell-Huey pontoon-equipped 
helicopters, shown in Fig. 11-2, which landed on the water at sampling stations 
ranging throughout the reservoir. Care was taken tlat the sampling crews con- 
ducted their surveys to coincide with the overpass of the Landsat satellites so as 
to Insure the validity of the data comparison. Each sample site was identified 
and marked on maps by sighting on prominent land features so that the exact 
location could be resampled on the rext Landsat flyover date. Samples were 
analyzed after acquisition at the Environmental Monitoring and Support Laboratory 
at Las Vegas and data recorded by the STORET system. However, due to weather 
conditions which restricted helicopter use and mechanical problems with the 
helicopters, sampling was severely restricted, as shown in TaLie I I — 1 . 
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Figure II-l. Map of Lake Mead and surrounding area 
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Figure II-2. An EPA Bell-Huey helicopter prepares to land on a lake to obtain water samples 
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Tabic 11-1. Landsat acquisition dates and corresponding 
ground truth sampling 





Sampled 

Date, 197b 

Frame No. 


by Kl’A 

June 1 1 

5419-18584 



June 12 

5420-17042 



June 20 

2515-17254 



June 2 1 

2518-17312 



June 29 

5437-18572 



June 30 

54 38-1 7030 



July H 

2533-17251 



July 9 

2534-17305 


X 

July 27 

2552-17302 



August 4 

547 3-18542 


X 

August 5 

5474-17000 


X 

August 1 j 

2598-17241 



August 22 

5491-16525 



August 2 1 

5492-1658 3 


X 

August 31 

2587-17234 



September 1 

2588-17292 



September 18 

2 GO 5-17231 



September 19 

2 606- 1 7 26 5 



September 27 

5527-16495 


X 

September 28 

5528-16553 


X 

The lack ot ground truth 

data lor all but 6 flyover dates 

limit 

ed accurate corro- 

lat Ion ol o lass i t teat ton 

results vs trophic Indicators to 

on 1 v 

6 Lind sat scenes. 

Of the six scenes, tour were cloud-free over the lake surface area while two were 

obscured by cloud cover. 

as Indicated In Table 11-2. 




Table 11-2. 

l.indsat 1 tame number, date, and approximate cloud 
cover tor sample acquisition periods 

Cloud 


Data, 197t> Frame No. Covet 


July 9 

2534-1 7305 

C 1 ea r 

August 4 

5473-16542 

C 1 ea r 

August 5 

1 

5 

C' ear 

August 2 3 

5492-16583 

Obscured 

September 27 

5527-16495 

C 1 ea r 

September 28 

5528-16553 

Obscured 
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4. 


Image Process In# Methodology 


The I PI- lias established a set of standard procedures to be applied to Landsat 
Imagery before any follow-on processing Is attempted. Figure 11-3 depicts the 
processing sequence used on l.ake Mead Images. The preprocessing procedures 
consist mainly of cosmetic und geometric corrections. In terms of cosmetic pro- 
cessing, Iwtndsat Images are corrected for MSS line dropouts, such as slipped or 
missing ! lnes and other obvious visual defects In the MSS Imagery. In addition, 
Landsat-1 MSS Images have been plagued by a striping problem which Is the con- 
sequence of an Imbalance in the sensor detectors. If this phenomenon occurs In 
a selected Image, the frame Is preprocessed to reduce the effect of the striping, 
also known as sixth-line banding. 

Geometric corrections are achieved through the routine application of the 
V 1 GAR program VEKTSLOG. Data contained Ir. the Landsat CCT's, as supplied by the 
EROS Data Center, are In a four spectral band interleaved format which is not 
compatible with the processing approaches established at the 1PL. The iPL soft- 
ware program VKRTSLOG unravels these Interleaved data, producing u separate Image 
for each spectral band. A history label Is also produced which supplies the user 
with pertinent information about the Image, such us latitude, longitude, scene 
identification and uny processing history. The program also applies geometric 
corrections to the data, including corrections for mirror scan velocity, panorama 
correction and resampling of the data to create an approximately 80-meter in- 
stantaneous field of view (1F0V). 


5. Registration of Extracted Lake Mead Subscenes 

Due to the large number of Landsat Images of l,ake Mead which were to be used 
in the study, the resources of the IPL's library of over 300 documented image 
processing programs were used to aid in the development of a more easily wielded 
data base. The decision was made to extract subscenes of the Lake Mead image area 
from each frame and to register these subscenes to one reference image so tliat all 
pixel coordinates would be identical from scene to scene. 

Frame 547i-lb542 (August 4), depicted in Fig. 11-4, was chosen to serve as 
the reference image to which all other landsat frames would be registered. The 
frame appeared to be free from cloud cover and the entire lake surface was con- 
tained within the scene. The program P1CREG was used to collect tiepoints which 
correlated the Images to be registered. Figure 1 1-5 schematically represents 
the operation of this program. P1CREG Is an Interactive VICAR applications pro- 
gram designed to aid in the registration of two images. PICREG will display two 
pictures on a split-screen video display, with each picture window 512 lines by 
320 samples. The user positions trackball-controlled cursors over identical 
features in each image and enters these tiepoints by keying in a predefined cliar- 
acter. Several modes are available to improve the input tiepoint position, and 
the user has the capability oi recursor ing any bad points or of restarting at any 
posit ion . 

The correlation mode in PICREG was used, which attempts optimization of a 
tiepoint position using two areas of linear dimension 2^. The default value N * 5 
was used; thus the size of the correlation region was 32 pixels. PICREG requires 
that the user must choose at least five tiepoints. For the purposes of this 
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Figure 11-3. A reproduction of a poster board used to describe the image process- 
ing methodology used at JPL for lake water quality monitoring 
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Figure II-4. Lind sat frame 3473-16542, August 4, 1976, of the Like Meal Study 
area. This image has been contrast-enhanced prior to falsi- 
product ion 
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MAP PICTURE 1 ONTO PICTURE 2 BASEO CN A GEOM DERIVED BY 
LOCATING COMMON FEATURES 

Figure 1 1 — 5 . A graphic Interpretation of how the various Lake Mead Imaged were 
registered to each other 


project, governed by the size of the subscenes, at least 24 tlepolnts were used, 
depending upon the amount of lake surface area which was contained In a particular 
subscene . 

Tlepolnts were located on geologic formations for the most part, as these 
formations were highly visible and relatively unchanging. The tlepolnts were 
located within close proximity to the lake shoreline, but not on the actual shore- 
line. This procedure alleviated confusion associated with fluctuating water 
levels and sliadows at the shoreline. The attempt was made to use the same geo- 
logic formations from frame to frame, if the formations appeared in the scene and 
were not obscured by cloud cover or shadow. 

After the tlepolnts were located and before the geometric transformation was 
performed, PICREG was allowed to reorganize the tlepolnts through the "fit" option. 
The fit parameter allows PICREG to fit a surface, where (x,y) are the reference 
picture and (x’,y') are the input picture. 
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to all of the tlepolnts and t* Interpolate new tlepolnts Iron this expression Into 
u 5 x 5 matrix within the area encompassed by t he originally specified tlepolnts 
on the reference picture. At this point, the user could a 1 ho check the accuracy 
of the "fit" and either trv a n«.w approach or accept the fit ax computed. 

If the fit wuh deemed acceptable, within l pixel accuracy, the geometric 
transtorm.it Ion wax t hen performed. Following the transformation, the accuracy 
of the registration was checked visually by <. resting •* difference p.< lure using 
the reference Image and the transformed Image. All 20 Utndsat scenes were reg- 
istered to the control scene using PICREO. This created a data base in which 
the coordinates for each sample site would be Identical for each laindsat image. 


b. Sample Site Location 

The NES handled the problem of accurately locating sample s i t vs on t fie lake's 
surface by sighting the location of the helicopter as samples were gathered. The 
position was located by sighting on prominent land or cultural features with the 
helicopter compass and then calculating its distance t rom the shore and/or 
features. locational Jata were then recorded in the field notes and later entered 
Into the STORET systeri along with the trophic indicator data. 

Color aerial transparencies of Lake Mead, taken at altitudes ranging from 
23,200 feet to 23,900 feet using a Zeiss camera equipped with a 152.4 mm lens, 
were supplied by JPL to EPA/LV. EPA/LV, using the locational data from STORET, 
marked the locution of each water sample station with a pinpoint and the cor- 
recpondlng STORET Identification number and returned the photographs to JPL. 

Using these as a guide, the sample sites were then located on the reference image, 
frame 5473-10342 (August 4), with the aid of a bendlx Datugrld Digitizer. Figure 
I l-b shows tlie reference frame with sample sites outlined. 

The mat* «.x size of pixels to he used for the spectral analysis was 5x5 for 
all but four of the 49 total sample sites. Sites VR04 , LV08, LV09 and IC01 were 
located in are is cf Lake Mead in which the shoreline configuration necessitated 
that the site boundaries be adjusted to fit within the lake's confines. Table 
11-3 gives the dimensions of the pixel matrices for these four sites. 

The Landaat imagery used by IPl. Is, in most cases, resampled to produce 
pixels measuring approximately 80 meters by 80 meters. More specifically, the 
resampling has produced square Land sat p';tels which represent an earth surface 
distance of 79.98 meters per pixel edge. This resampled pixel represents an 
area of b,39b.8 meters- (b8,854.b feet-) or 0.b39b hectares (1.58 acres). There- 
fore, the selected 5x5 pixel matrix represents 15.99 hectares or 39.54 acres on 
the lake's surface. Figure 1 1 — 7 illustrates both one pixel and the 5x3 pixel 
matrix in relation to a standard U.S. one-mi ie-square section. 
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Figure 11-6. Lake Mead with the water sample site locations marked by small 
white squares 
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Table 11-3. Sample sites not conforming to the standard 5x5 matrix size 


Sample site 


Hilt t lx size 


VR04 

5 

X 

2 

LV08 

2 

X 

2 

1.V09 

1 

X 

5 

I an 

5 

X 

2 


Thu original Intent of the project was to classify Uike Mead solely using 
supervised classification routines which were developed for the previous lake 
tasks (Refs. 11-1,3,4). One limitation of this approach Is tliat these programs 
depend heavily upon statistics gathered from sample sites. Man? times sampling 
datu for a lake are either out of date, simply too expensive to gather or not 
available. This necessitated the development of unsupervised routines which did 
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Figure 11—7. A graphic representation showing the area of Landsat pixel(s) 
in relation to a standard U.S. section 


3b 


ORIGINAL PAGE IS 
OF POOR QUALITY 



not rely on training sites. During the course of the project, the 1PL h» d become 
Increasingly involved In the development of unsupervised clustering routines. An 
unsupervised program, USTATS, wub applied to the August 4 scene 5473-16542. The 
results from the unsupervised routine were encouraging for this particular image, 
in terms of a comparison with the sample site measurements provided by EPA/LV. 

The unsupervised program had attractive benefits which enhanced its continued use 
by the project. lor example, the routine did not require the somewhat timely 
gathering of MSS statistics and ranking of classes associated with the supervised 
routines used previously. A lest of the supervised classification routine yielded 
classification maps which indicated confusion in class discrimination as indi- 
cated in Fig. ll-H. It was felt that if the unsupervised classifier could prove 
to be as accurate as or more accurate than the supervised techniques, more em- 
phasis would be placed upon its development and application to the Lake Mead 
study. 


7. St. tistlcal Lathering from Sample Sites 

Although interest in unsupervised classification Increased, it was decided 
to continue with procedures developed in earlier projects until the unsupervised 
technique could be refined and further investigated. To this end, sample sites 
which had been located and identified on the reference image were then located on 
.ill of the registered images through the application of the interactive program 
STATS2 . 

Follow-on supervised mult lspectral classification programs developed and 
used previously at IPL required the creation of MSS statistical data sets. The 
data sets contained pixel values for each "training" (or sample) site to be used 
in an MSS classification routine. The introduction of an interactive image 
processing system at IPL enabled programmers associated with this project to 
develop a fast routine for accurately locating training sites, computing the 
associated statistics and storing the information in a usable form. The program, 
STATS2, detailed in Section I-C, enabled the user to display a subject image of 
variable size on a video screen. The image was input in MSS format (interleaves 
and contained up to four spectral channels. The image could be displayed in any 
one of the input channels and was viewed through a movable "window" of 512 lines 
by 640 elements. The program also offered a magnification option which enabled 
the user to accurately locate a training area and vary its size to remain within 
the boundaries of the subject area, such as a lake shoreline. The user could 
select training sites by either keying in predetermined coordinates or by locating 
a site visually and defining it through manipulation of a trackball controlled 
cursor . 

Coordinates for this project were keyed in for the preregistered Landsat 
images and assigned a 5 x 5 matrix pixel density. The matrix size was charged 
from 5x5 only if the sample site overlapped the lake boundaries, as in the case 
of sites located in a few narrow areas. After a site was located and drawn, a 
command would be given to the program to display the statistical information for 
the site. Statistics were returned for all input channels in the form of x and 
sigma and the size of the pixel matrix was also displayed. A unique name was 
assigned by the user to the site at the time it was saved, which served as iden- 
tification in future reference. The NES STORF.T number was used to identify I.ake 
Mead water sample sites. Once a site was saved, the user returned to the proces- 
sor to continue with the selection of the next site. 
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Figure II-8. A supervised classification image of Lake Mead using a principal 
component transformation of 10 trophic indicators as a basis to 
effect classification 
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Although eventually the sample Hite Ntatlatlca were not heavily relied upon 
for claaalf leal Ion, the statistics did provide a moat Important factor In the 
discovery and resolution of data discrepancies which plagued the project. 

8. Data Discrepancies 

As the MSS statistical data sets were being created by STATS2, comparisons 
were made with day-to-day spectral signatures to check the reliability of the MSS 
data. When spectral signatures from sample sites In the June 11 frame 5419-16584 
were compared with signatures from the following day's coverage (5420-17042), the 
presence of serious Inconsistencies became obvious, as Illustrated In Table 3. 

As the operat Ion of STATS2 continued on subsequent images, further comparisons 
were made and t lie problem was found to be occurring to some degree in all consecu- 
tive day coverage. In the case of a comparison between s Ignat ure« on August 22 
frame 5491-16525 and August 23 frame 5492-16583, the discrepancy ranged up to 
10 DN, as seen In Table 11-4. 

Several possible causes for the Inconsistencies were considered and rejected. 
(Examination of the EPA water quality measurements, for conflicting sample sites on 
consecutive days, did not support a theory that the water quality had changed 
enough to seriously affect the MSS spectral response. Atmospheric changes were 
also considered. The effects of clouds and atmospheric haze on the spectral 
response of l^indsat sensors has been documented (Kefs. 11-7,8). However, a 
comparison of Images which contained cloud cover with cloudless scenes did not 
entirely support the atmospheric disturbance theory. Further, as suggested by 
Kverctt, Leon hurt and l.epley in their Investigation (Ref. 11 — 1), the remoteness 
of the Lake Mead area serves to minimize the probability of atmospheric dis- 
turbances related at least to human activities. This was precisely a governing 
reason for selecting the area, both In theii study and In this project. Had the 
data been acquired by both l.andsnt-1 and Landsat-2, a theory regarding the possibl 
difference between the spectral response ol the two spacecraft sensors would have 
been investigated. However, ail consecutive day coverage was provided by the 14t 
overlap In scan paths from l«tndsnt-2. 

Having virtually ruled out all suggested considerations, It was concluded 
tliat 1) the discrepancies could have been caused by variations in the sensor 
response which were not detected and therefore not corrected in the ground cali- 
bration procedure or that 2) atmospheric changes had occurred which were unknown 
and undetectable. 


9. Resolution of Data Discrepancies 

Regardless of the cause, the discrepancies had to be corrected so that the 
project could continue. In deriving a solution to the problem, two cons lderat Ions 
were of primary concern: 1) the proposed correction should not degrade the data 

quality; 2) due to time considerations and cost constraints, the correction hau 
to be efficient and simple enough not to require large amounts of comput-r 
processing time. 
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Table 11-4. 


examples of discrepancies 


Date, 197b 

Site 

Band channel 

DN 

June 11 

0A02 

4 

45.1 

June 12 

0A02 

4 

40.1 

June 29 

BB05 

4 

50.2 

June 30 

BB05 

4 

44.9 

June 29 

0A02 

4 

57.8 

June 30 

0A02 

4 

53.7 

August 22 

BB05 

4 

46.7 

August 23 

BB05 

4 

36.6 


Three ideas were proposed and judged on the basis of these two concerns. 

The first suggestion, to attempt to recalibrate the data using the calibration 
wedge which is provided with the landsat imagery, did not sati? y criterion #1. 
Previous attempts to use the calibration wedge at iPL for similar purposes had 
proven Ineffective. The calibration wedge provided with the landsat CCT's con- 
tains only six sample points, of which two are typically saturated. The remain- 
ing four points do not provide enough information to properly recalibrate the 
data. It was felt that further data quality degradation was highly probable with 
this method. 

The generation of an atmospheric model, wnich could account for some of the 
error, was also considered. Work has been performed in this area by other in- 
stitutions with a degree of success (Ref. 1 1-9). However, it was felt that this 
approach did not satisfy the time constraints of criterion #2. The procedure 
would have to be applied to 20 separate landsat scenes. The computer processing 
Lime would have been costly in both hours and money. 

The third alternative was to devise a normalization procedure which could 
operate on the MSS (interleaved) data in its present, form, correcting all four 
spectral bands in a single operation. This hopefully would satisfy the time 
constraints. It was thought that by selecting a "control area" from one landsat 
scene which was of relatively unchanging spectral response, it could be used to 
essentially normalize or conform the identical area in the remaining landsat 
Images. Similarly, all pixels comprising each Landsat frame would then be 
normalized, based upon the differences between the standard frame's control area 
and the same area in the subsequent frames. By constraining the spectral signa- 
tures to be constant in the control area, it was believed that a correction would 
be made for both the sensing system and any atmospheric effects simultaneously. 
Initial tests Indicated tliat this procedure should be selected due to the ease of 
implementation and the results it could produce. 

The control area, as seen in Fig. 1 1-9, was selected in the Virgin Basin area 
of Lake Mead. A comparison of several test sites ranging throughout the lake indi- 
cated the Virgin Basin's spectral response to be the most stable in all the Landsat 
scenes. The spectral signature for the Virgin Basin Control area was calculated 
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Figure 11-9 


A false color image of Lake Head showing Che upproxliaate location In 
the Virgin Baaln where MSS data for 13 daya were ccjpared 
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for each of the Landsat scenes. The August 4 image 5473-16542 was selected as the 
standard control area to which all other scenes would be normalized. The August 4 
scene was cloud-free and did not exhibit any haze. The DN differences between the 
standard control area and a given t rame were then added to or subtracted from each 
picture element in the given frame. This procedure proved to be effective, inex- 
pensive and easily implemented into the project plan. 

The effect of the normalisation procedure performed on the Lake Mead scenes 
wus analyzed in two ways: 1) the statistical effect on the data was analyzed, and 

2) the effect on c lar sif icat ion performance was examined. 

In the first case, the standard deviation of selected sample sites from all 
the scenes was computed for each spectral channel using STATS2 . The standard 
deviations dropped significantly after the normalization was performed. As an 
example, the standard deviation for the ilverton Arm sample site (0A02) dropped 
from 9.8 to 1.7 in band 4. Line plots (Fig. 11-10) of the sample sitn DN's vs. 
scene (date) illustrate the effect the normalization had on stabilizing DN 
variation from scene to scene. Figure 11-10 depicts line plots comparing Landsat 
data before and after normalization was applied. 

To test the effect on classification performance, an unsupervised classifica- 
tion was performed on both June 11 (5419-16584) and June 12 (5420-17042) scenes. 
Using an unsupervised routine, clustering of the water signatures was performed 
on the June 12 frame before normalization was applied. The statistics for these 
clusters were applied to both the June 11 and the June 12 data. The same pro- 
cedure was then followed for the normalized output of the same frames. The 
results Indicated that the normalization had caused the June 11 scene to classify 
more closely to the June 12 scene, as shown in Table 1 1 — 5 . 


10. Final Data Preparation for Classification 

The data discrepancies having been resolved, the project moved into the final 
pliase of data preparation prior to classification. At this point, each image sub- 
scene contained spectral information for the water surface area of Lake Mead as 
well as the surrounding land features. Since only the water data was of Interest, 
the land surface area information was "removed" from the image before classifica- 
tion. This allowed the clustering routine to operate most efficiently, by avoiding 
confusion and unnecessary computer time in water/land discrimination. The resources 
of the LAKELOC program were used to isolate Lake Mead from the surrounding land 
mass. Figure 11-11 depicts a Lake Mead subscene in the final binary mask form as 
output from LAKELOC. 


11. Lake Mead Classification Summary 

The thematic maps illustrated in this section were all derived from the appli- 
cation of USTATS. USTATS is an unsupervised routine which accumulates statistical 
information about a scene using local clustering. Results from unsupervised 
classifications using USTATS were more satisfactory than those obtained from super- 
vised techniques used in the past. The algor ithm implemented in USTATS was orig- 
inally developed at Pennsylvania State University under the name CLUS, and has 
undergone extensive modification to run in the VICAR system. After USTATS lias 


45 


PRECEDING PAGE BLANK NOT FTLMEB 




A. 

W 




8. 

< n 



O 

D 

< w> 


o 

D 

< V 


D 

— i o 


Z 

3 


8 


z 

3 & 


z 

o 


8 


z 

D tv 


z 

3 - 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Figure 11-10. A plot showing raw and normalized MSS means and standard deviations for the green, red 
and IR1 channels. Data shown are from the center of the Virgin Basin for the dates 
indicated 



Table 11-5. Effects of normalization on classification performance 


Percentage distribution by class 

Class 

No. 

June 12 

June 11 
(Original) 

June 11 
(Normal lzed) 

1 

65.2 

16.2 

55.0 

2 

16.7 

37.3 

16.2 

3 

1.5 

14.0 

5.1 

4 

7.5 

4.8* 

0.6* 

5 

6.9 

17.6 

10.6 

6 

1.1 

3.3 

2.6 

7 

1.1 

0.5 

3.6 

Classification affected by clouds. 


completed the cluster process, it passes statistics to the VICAR program BAYES, a 
Bayesian maximum likelihood classifier. baves performs Lhe actual classification 
using a non-Euc 1 ldean distance function to decide which cluster signature is most 
similar to that of a given pixel. Figures 11-12, 13, 14, and 15 are thematic 
photomaps of Landsat scenes for August 4, August 5, July 9 and September 27 re- 
spectively. These four scenes are the only ones which have supporting ground 
truth measurements which can be used to compare classification r. curacy. 

Visual examination of these naps indicated they provide llmnologlstn with an 
excellent synoptic view of I.ake Mead. In addition, the areal extent of river plumes, 
as seen in the Overton Arm and Pearce Ferry areas, was successfully mapped. Spec- 
tral classes were correlated most effectively with Secchi depth measurements, one 
percent light level and nephelometric turbidity. These indicators were expected to 
correlate well due to their direct correspondence with visual phenomena produced by 
turbidity and thus affect the spectral response. Chlorophyll content was also fairly 
successfully correlated with spectral classes. Complete analysis of the classifi- 
cation results is presented within the following section. 

Two deficiencies evidenced by previous lake classification projects, conducted 
by the 1PL, were encountered again by this pioject. Visual correlations between 
spectral classes with total phosphorous and total organic nitrogen measurements were 
difficult to construct. The sixth-line banding problem which has plagued Landsat 
data, and the data users, continued to degrade the resulting classification pro- 
ducts, despite attempts to minimize its effect. 




Figure 1 1 — 1 1 . A binary mask of Lake Mead 


Despite these deficiencies, l.indsat demonstrated great potential for use In 
the monitoring of changes In lake Mead. This is evidenced by an examination of 
Table ll-b, which correlated the results from the unsupervised classification of 
four l-andsat scenes with corresponding water quality measurements made f . om sur- 
face samples on the same Landsat flyover dates. 
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Figure 11-12. Unsupervised classification of Lake Mead for August 4, 1976 




Figure 11-13. Unsupervised classification of Lake Mead for August 5, 1976 
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Figure 11-14. Unsupervised classification of Lake Mead for July 9, 1976 
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Figure 11-15. Unsupervised classification of Lake I. for September 27, 1976 
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12. Analysis of Classification Results 


Reference la made In the following text to limnological measurements and 
conditions which may be unfamiliar. Briefly, the productivity of an ecosystem, 
such as a lake, connotes a quality whereby living substance is manufactured tnrough 
interactions of community and environment (Ref. 11-6). The productivity or 
fertility A a lake depends on nutrients received from regional drainage, on the 
geological age, and on the depth of the water body (Ref. 11-10). '1'u‘bldlty is 
the degree of opaqueness produced by suspended particulate matter (Ref. 11-6). 
Secchi Disc transparency (SDT), one percent light level (IX level), nephelometric 
turbidity (NT) and chlorophy 1 1-a (CHLA) are measurements which are indicators of 
productivity and turbidity. 

The statistics, illustrated in Table 11-6, were derived from the correlation 
of the thematic maps with contact-sensed data in the following manner. An image 
polygon output from STATS2, containing the 5x5 pixel sample site matrices, was 
overlayed on each thematic map. Each sample site matrix was then examined for 
its homogeneity as evidenced by the color or colors represented within it. Sites 
located within highly heterogeneous matrices were dropped from the analysis. Al- 
though it is recognized that a pixel measures HO x 80 meters as processed by the 
1PL and tliat a sampling site is readily contained within a pixel, a degree of 
uncertainty exists as to which pixel within the matrix is the one containing the 
sampling point. Thus it was decided to drop all highly heterogeneous sites. All 
remaining sites which were identified as falling within a particular color-coded 
spectral class were then pooled and descriptive statistics generated for earn 
ground truth parameter illustrated. The mean measurement is listed first, followed 
by the range of ground truth measurements for that class. The last statistic, N, 
equals the number of sampling sites used to calculate the mean. The statistics 
generated were then used to characterize the spectial class in terms of water 
quality. 

Spectral classes from the thematic maps can be pooled into three general 
categories. The first category, PRODUCTIVITY, is represented in the thematic maps 
by dark blue, light blue and green. Areas color-coded dark blue can be described 
as ol igotrophic , as evidenced by the SDT, low CHLA and deep 1% Lev. measurements. 

e light-blue-coded class can be considered to be a borderline between uligo- 
trophic and mesotrophic, based upon supportive contact-sensed measurements. The 
third productivity-related class does not appear in all thematic maps. It is a 
relatively small class with little supportive ground truth. However, based upon 
SDT and 1% Lev. measurements and the geographic locations of the class, it appears 
to be more productive than either the dark blue or light blue class. 

The second category Iras been designated MIXED EFFECTS and is represented by 
a gray class occurring in some thematic maps and a black class appearing in all 
the images. It is suspected that the gray class does not relate to water quality 
but is evidence of residual sixth-line banding which was not successfully removed. 
The black coded class lias no corresponding water quality measurements. However, 
an examination of the lake in terms of where the black pixels occur, coupled with 
a knowledge of the lake, strongly suggests that the class consists of the inter- 
mingled effects of nearshore areas. These areas typically exhibit effects caused 
by visible lake bottom, submerged plant cover and shadows. 
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The color coded areas of yellow, brown and orange represent the third cate- 
gory, TURBIDITY. The yellow coded areas are characterized by a lowering SDT and 
Increased NT. The brown coded areas can be readily identified as river water 
which has entered the lake and Is still In the process of dropping Its suspended 
sediments. The most turbid class of water found In Lake Mead Is represented In 
the orange coded areas. These areas consist of the heavily sediment laden 
plume of the Colorado River where It enters Fierce Basin and a small portion of 
the Muddy River branch of the Overton Arm. These areas signal the influx of a 
substantial quantity of inorganic materials into the lake which effectively reduce 
the zone In which photosynthesis occurs, thereby limiting productivity. 


B. THE LAKE TAHOE STUDY 

In order to investigate the possibilities of a comprehensive data base sys- 
tem which took into account the contributing factors as well as the effects of 
lake eutrophication, the IPL chose to construct an IBIS data base for the Lake 
Tahoe Basin. Due to increasing environmental pressures and rampant land develop- 
ment, this region has been the subject of intensive study by the U.S. Department 
of Agriculture and the Tahoe Regional Planning Agency (TRPA) . In 1977 TRPA made 
available a volume entitled. Water Quality Prob lems an d Management Program , which 
detailed present knowledge of the lake Tahoe Basin Watershed. Using data compiled 
in this volume as a source, analysts at IPL proceeded to integrate these data 
with Landsat and digital i. rain data produced by the Defense Mapping Agency 
(DMA). The IBIS data base method allowed cross correlation of Landsat imagery 
and topographic data with USDA environmental data relating to such parameters 
as drainage basin acreage, surface runoff and terrain configuration. Data base 
construction led to the production of integrated images created by "overlaying" 
conventional maps and digital maps on Landsat. Preliminary work was also per- 
formed which attempted to depict the terrain configuration for each drainage 
basin based on DMA topographic data. 


1. The Lake Tahoe Basin 

The Tahoe Basin occupies over 500 square miles (1295 km^) situated in a 
graben straddling the boundary between California and Nevada. Figure 11-16 re- 
produces a map of the Lake Tahoe Basin and its environs. Lake Tahoe contains 126 
million acre feet (155.4 km^) of water in a 190-square-mile (492.1 km^) surface 
area receiving inflow from 63 tributaries with only one outlet at the Truckee 
River. The Tahoe Basin has traditionally attracted recreational activity due 
to its clear deep water and pine-forested shorelands coupled with its proximity 
to major metropolitan areas in northern California. Since the 1950's the basin 
lias experienced escalating demands for land development at the expense of the 
natural watershed. Discharge of sediments to the lake has greatly increased due 
to accelerated human interference, and alterations to the natural drainage patterns 
are evident in some areas. According to the Water Quality Plan published by the 
Tahoe Regional Planning Agency, although most development lias occurred in areas of 
mixed-pine and pine-fir vegetation types, the greatest impact of urbanization upon 
vegetation and hence the watershed environment has occurred on several marshes 
and meadows within the Basin. 
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Figure 11-16. Map showing Che geographic location of Lake Tahoe 


The problems which the Lake Tahoe b'sin is presently confronting are certain- 
ly not unique to this area. The consequences of man's alteration of the natural 
environment are symptomatic of the pressures which can be brought to bear by 
the activities of an increasingly mobile populace seeking new recreation and 
living space. 

2. Construction of the Lake Tahoe Basin Data Base 

a. Land sat Data. IPL chose to construct the Lake Tahoe Basin data base 
from data sources which were readily available and familiar to its analysts. 
Landsat imagery, which lias formed the basis for the majority of IPL's water qual- 
ity assessment efforts, was chosen as the planlmetric base to which all other data 
would be registered. A subsection of a landsat 2 scene from August 27, 1976, 
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was extracted and t ransfor med Into a lambert conformal conic projection to serve 
as the planlmetrlc base as depicted in Fig. 11-17. To this scene two other Land- 
sat scenes were registered, one from Landsat 2, July 21, 1978, and one from Land- 
sat 3, July 3U, 1978. Summer scenes were chosen to reduce the chance of snow 
covet Interference. Registration was achieved through the aid of a piecewise 
surface fitting algoritlun. The surface transformation was defined through a 
series of tiepoints selected during the use of an Interactive spatial pattern 
recognition routine (PICREG) in which the analyst selected common geographical 
features (Refs. 11-11, 12). (Details of the PICREG program can be found in the 
technical application section on lake Mead). After all Images were registered, 
color reconstruction was also performed on each scene using a color enhancement 
technique which approximates a gauss lan distribution using the principal com- 
ponents of each landsat image separately (Ref. 11-13). This type of enliancement 
was designed to attempt even distribution for all regions within an image to 
avoid the saturating effects of conventional linear enliancement s. 


b. DMA Digital Topogra p hic Data . Terrain data produced by the Defense 
Mapping Agency (DMA) was next integrated into the data base. Digital terrain 
tapes were acquired for the Like Tahoe region from the National Cartographic 
Information Center. The tapes were prepared from U.S. Geological Survey 1:250,000 
scale topographic quadrangle map series. Map contour lines falling within a 
1-degree block of latitude and longitude were digitized and a matrix of elevations 
generated with one elevation tor every 0.01 inch on the map (200 feet/60.96 
meters on the ground). The terrain data were reformatted for use in IPL's 
VICAR operating system. In reformat t ing, lialfword integer elevation values were 
converted to single byte integers and scaled to the terrain variation within the 
area (Ref. 11-14). 

After reformatting, the terrain data sets were rotated 90 degrees counter- 
clockwise to orient the north-to-south values vertically in the sample direction, 
placing north at the top of the images. This compensated for data format as 
produced by the National Cartographic Institute. For the Lake Tahoe Basin, four 
separate terrain quadrangles were required to construct a complete image compris- 
ing the Tahoe Basin area. This necessitated the mosaicking of the four quadrangles 
before final registration to the data base. Figure 11-18 reproduces the four 
digital terrain quadrangles required to completely encompass the Lake Tahoe Basin 
and the final mosaic image. The relief-like effects portrayed in the Images were 
produced by digitally shading between contour intervals. 

Registration of the terrain mosaic with the Landsat planlmetrlc base was 
achieved through the application of a resampling algorithm which applied a two- 
dimensional correction grid derived from selected control points (Refs. I I —1 1 , 12). 
Tiepoint selection was achieved interactively (PICRF.G) relating surface features 
through comparison of the Landsat planlmetrlc base image and a relief-like version 
of the terrain image. 


c. Grap hical Data. Graphical data for the lake Tahoe Basin data base were 
acquired from conventional maps produced for the Tahoe Regional Planning Agency 
by the U.S. Department of Agriculture Forest Service. These maps presented hydro- 
logic and climatologic data pertaining to the lake Tahoe Basin. Before integration 
into the IBIS data base, graphical data must be transformed into image foimat. 
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Figure 11-17. False color composite image of Lake Tahoe from Landsat frame 5496 
17202, August 27, 1976. 
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4-1 : 250,000 

DIGITAL TERRAIN QUADS 

Figure 11-18. Four digital terrain quadrangles were mosaicked to produce a ter- 
rain map which contained the Lake Tahoe Basin 
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ThU was achieved by first digitizing the graphical data onto magnetic tape. 
Next, a least squares affine transformation was applied which created a general 
surface fit. Corresponding line and sample coordinates, In terms of the landsat 
planlmetrlc base, were calculated from latitude and longitude coordinates taken 
from the graphical map data for geographical features common to both data sets. 
Registration was once again achieved through the Interactive (P1CREC) selection 
of control points which linked geographical features from the graphical d«Ca file 
and the planlmetrlc data base. In Fig. 11-19, the transformation of a convention- 
al map to a graphical data plane Is Illustrated. Step A represents the digitizing 
procedure In which the drainage basin map Is converted to graphical Image format. 


d. Georefere nce Plane. A georeference image plane which provided an inter- 
face between all dataplanes for the Lake Tahoe Basin data base was created from 
the drainage basin map provided by the Tahoe Regional Planning Agency. As seen 
In Fig. 11-19, this map identified the location and outlined the boundaries for 
the 63 drainage basins Identified as providing watershed inflow to Lake Tahoe. 

The map was digitized and subsequently registered to the planlmetrlc base. Each 
separate drainage basin was then assigned a unique data number for all picture 
elements comprising that basin. Step B, In Fig. 11-19, represents the process of 
encoding each drainage basin, a procedure known as "painting." The final product, 
the georeference plane, is reproduced at the right in Fig. 11-19. All tabular 
data corresponding to the drainage basins were then entered into the data base by 
referencing the unique number assigned to each basin. The georeference plane Is 
an integral element within the data base, for it is through this plane that all 
tabular and most Image data are Interfaced. Figure 11-20 is a conceptualization 
of the georeference plane for the Talioe region. The magnified portion of the 
figure Indicates the coded picture elements which comprise each basin. In this 
figure a most basic application is illustrated, in which the picture elements for 
each basin are summed and transferred to an interface file in which the sums are 
stored according to basin code number. This Interface file can be accessed at a 
later date to produce statistical output such as acreage estimates for each 
drainage basin. 

Figure 11-21 illustrates the Lake Tahoe Basin data base as completed, forming 
a series of registered data planes. Once all data planes are in registration and 
at least one georeference plane is registered, the data base can be used to develop 
additional products through manipulation and integration of the planes. 


3. Application of the IBIS Concept to the lake Tahoe Basin 

a. Integrated Ima ges . Manipulation and integration of data planes compris- 
ing the data base is achieved through the implementation of VICAR and IBIS programs 
developed aL 1PL. Standard products output from IBIS Integration consist of over- 
lay images. Figure 11-22 reproduces a precipitation map produced by the Tahoe 
Regional Planning Agency. Using IBIS resources, this map was electronically digi- 
tized and registered to the planlmetrlc base image. By adding the digitized map 
to the color version of the Landsat scene for August 27, 1976, an Integrated image 
is created. As seen in Fig. 11-23, this type of product can represent more clearly 
the relationship between mean precipitation and the Tahoe environment as depicted 
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Figure 11-19. iransformat ion of a conventional map to a graphical data plane and georeference plane 




GEOREFERENCE PLANE 
ORAINAGE BASINS 


Figure 11-20. 


Conceptual drawing of the georeference plane and tabular Interface 
file 
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A conceptual diagram illustrating the development of the Lake Tahoe 
data base from digital terrain data, conventional data and landsat 
imagery 
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Figure 11-22. Conventional map showing mean annual precip- 
itation contours over the Lake Tahoe Basin 
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Figure 11-23. Image showing the registration of precipitation data from 
a conventional map with l<andsat Imagery. 
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by Utndsat. In future applications the mean mnual precipitation laohyeta could be 
converted to a continuoua surface image. This process la almllar to the conversion 
of USGS topographic maps to terrain relief images as depicted previously In Fig. 
11-19. The precipitation surface image could then be cross-tabulated with Indi- 
vidual drainage basins, using tlie geort Terence plane, to produce estimates of 
runoff coefficients. 

Simllarily, another Integrated image was produced by combining digital ter- 
rain data with Iamdsat imagery. The Defense Mapping Agency (DMA) mosaic was 
transformed to display 600-ft. contour Intervals as illustrated in Fig. 11-24. 

The digital terrain image was then Integrated with the August 27, 1976, Landsat 
scene to produce the image shown in Fig. 11-25. Integrated imagery of this type 
can be used to facilitate interpretation of terrain elevation in terms of a I.and- 
sat scene. DMA data were also used to calculate estimates of azimuth and slop* 
for each drainage basin. Details of this procedure are described in the 
following section. 

The georeference plane was integrated with each Landsat scene to depict the 
boundaries of each drainage basin mapped by the TRPA. Figure 11-26 shows each of 
three Landsat scenes contained in the data base as it appears integrated with the 
drainage basin map. While imagery of this type is useful for cursory Interpreta- 
tion and display purposes, the ability to extract additional information from 
data plane manipulation is the key to the application of the IBIS concept. 

b. Data Plane Ma nipula tion. Manipulation of the georeference plane, in 
this case the drainage basin map, provides the opportunity for the analyst to 
more closely examine the unique charac ter 1st ic3 of each drainage basin as revealed 
by the Landsat sensors. By referencing each uniquely coded drainage basin in the 
georeference plane separately, the analyst can extract and display any or all 
basins as individual images. This application is illustrated in Fig. 11-27, in 
which the georeference plane was used to extract the Trout Creek Basin and the 
Upper Truckee Basin, which are displayed in color as separate images. If a multi- 
spectral classification is performed on the data, information can be similarly 
reproduced in a basin-by-basin format. Tabulations of land cover type can then 
be expressed for all ba3lns collective' y or each basin individually. 

The georeference plane can also be interfaced with tabular files to produce 
a statistical output such as that reproduced in Table II-7. In this table, the 
area of each uniquely coded drainage basin comprising the georeference plane was 
determined by summing picture elements and converting the sum to total acres, 
producing a Landsat acreage estimate. These data were then cross-referenced to 
tabular files which contain acreage estimates produced by the TRPA. All estimates 
were then ranked and output in the form of a computer listing which provided a 
comparison of Landsat vs conventionally acquired acreage estimates. 

The inclusion of DMA digital terrain data provided an opportunity to attempt 
integration of an Important element of ancillary, ground-based information with 
remote sensing data. The DMA developed these data by Interpolating existing USGS 
1:250,000 scale topographic mrps to produce ultrafine mesh digitized latitude, 
longitude and elevation contour data. As described earlier, four DMA quadrangles 
were digitally mosaicked to form a unified, continuous surface elevation image. 
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Figure 11-24. Digital elevation contours produced by processing 
digital terrain data 


7a 


ORIGINAL PAGE 

BLAGK AND WHITF P MOTHER A PH 





GRN P7AUG76 


N38-39U1P0-1? 
BO.OVI'PXL 
G£»10 - DS4 - I 


SUM LL 4 3 A2114 P*C1£1 SC«J[ 


PCGSTR - re INI ♦IMP 


ftVSxTAHRSMX 


Figure 11-25. Image showing the registration 
of elevation contours with 
Land sat l.aagery. 
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LANOSAV SCENES CONTAINED 
Wl HIN THE IBIS DATA BASE EACH 
INTEGRATED WITH THE GRAPHICAL 
DATA PIANE VERSION OF THt 
DRAINAGE BASIN MAP 


Figure 11-26. Integration of data planes produces overlay image 
c reaving new visual tools which combine data from 
various sources and for different dates. 
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Figure 11-27. Interfacing the georeference plane with I.andsat imagery allows ac- 
cess to each basin separately and the generation of tabular reports: 

(a) - Trout Creek Basin 

(b) • Upper Truckee Basin 

This Image was then geometrically registered with the Lindsat Imagery to permit 
future cross-tabulations of elevation Information with other Information data 
planes within the data base. In an effort to quantify nnd process the digital 
terrain data, the elevation Image was processed with VICAR software to produce a 
component representing slope magnitude. To develop this component it was necessary 
to compute the vector crossproduct between the horizontal (east-west) Image 
elements and the vertical (north-ciuth) picture elements. This vector product 
then provided an estimate of the slope between adjacent east-west and north-south 
elevation Image elements. The outputs of this process were then coded to reflect 
slopes between 0® (no slope) and 90® (vertical slope). The angles In degrees were 
further coded for image output by rescaling 0° to be equal to 0 digital number (ON) 
or black and 90® to be equal to 255 DN or white. Figure 11-28 illustrates an 
application of these concepts. First, the slope magnitude was cross-tabulated with 
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Table 11-7. 


Lake Tahoe drainage oasina acreage tabulation 
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GEOREFERENCE 

PLANE 


DIGITAL TERRAIN 
GRADIENT MAGNITUDE 
IMAGE 


DIGITAL 

TERRAIN 

INTERFACE 






GRAOIENT MAGNITUDE IMAGES BY DRAINAGE BASIN 

Figure 11-28. Using the drainage basin map as the georeference image, the 
slope and gradient magnitude for each drainage basin may be 
displayed and tabular information produced 
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the da Inane basin georef erence plane. Each alope magnitude Image element associ- 
ated with the Individual drainage* basins was then extracted. Further cross- 
tabulation permitted the computation of u mean slope estimate by averaging the 
slope magnitude Image elements tor each drainage basin. 

Table l 1-8 reproduces an IBIS table generated from the digital terrain 
Imagery which lists azimuth and slope for each drainage basin. F.ventually, It Is 
hoped to develop a model for drainage basin terrain based on digital Imagery. At 
present, however, difficulties have been encountered with discrepancies between 
Individual DMA map quadrangles, especially along map edges, which have precluded 
the development of accurate models. Other types of digital terrain Imagery are 
being Investigated with the hope of Integrating more reliable data Into the lake 
Tahoe Basin data base. 

Table 11-8. lake Tahoe drainage basins azimuth (aspect) and 
slope magnitude statistical summary report. 
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A data bane, no matter how easily constructed and manipulated, la only as re- 
liable as the elements of which It Is made up. Therefore, as new data are acquired 
It Is essential that these be tested and Implemented Into the data base to Insure 
Its viability. Although existing digital terrain Imagery luis proved difficult to 
adapt, Landsat data continues to work well within the data base concept. However, 
in terms of water quality, the subject lake and watershed area must be of such 
dimensions as to accommodate the resolution of the Landsat sensors. 

Accurate correlation of changes In Lake Tahoe with various data elements con- 
tained within the data base will require study < ver several years. It Is hoped 
tliat such a system over time can be used to monitor and evaluate causes for 
changes which affect the lacustrine environment. This will require the development 
of precipitation modelling, surface runoff models and classification of drainage 
basin cover types. These elements must in lum be integrated and evaluated for 
accuracy before the system can be considered usable. Such a system Is feasible 
given the continued Improvement of the remote sensing tools used to construct the 
data base and the data Integrated Into it. 

C. STUDY OF SELECTED LAKES IN MONTANA, UTAH AND MICHIGAN 

Four ureas were chosen for study as a final demonstration of lake classifica- 
tion technology developed at JPI.. In Montana, Flathead l^ke, located in the 
northwest area of Lite slate, and several lakes within the Kootenai National Fores', 
nearby were selected. The Flaming Gorge Reservoir astride the northern border of 
Utah and southern Wyoming was also chosen. The Illinois State EPA also expressed 
a desire to be included in the final demonstration, agreeing to coordinate sampling 
efforts, on the same lakes studied previously (Ref. 1*3), with landsat flyover 
dates during the summer of 1978. The U.S. Forest Service, Ottawa National Forest, 
also expressed an Interest in surveying several lakes within their area in the 
northern part of Michigan. Therefore, In April and May of 1978 the IPL attempted 
to coordinate ground sampling of selected lakes In these areas with the USEPA 
Western Energy Office in Denver, Colorado (Kootenai National Forest and Flathead 
Like and Flaming Gorge Reservoir), Illinois State EPA in Springfield, Illinois 
(selected Illinois lakes) and the U.S. Forest Service, Ottawa National Forest, in 
Ironwood, Michigan (selected Michigan lakes). The intent was to once more test 
and refine the techniques for lake isolation and classification which had been 
under development since 1974. 


1. Data Acquisition 

One of the most significant problems associated with any study using Landsat 
data is with data acquisition. The problem is basically twofold, with each part 
independent of the other. l,andsat scene availability is always dependent upon 
favorable weather conditions; cloud cover can virtually rule out the use of any 
number of Landsat scenes. The other problem, assuming weather conditions to be 
fair, is scene availability from the data source, such as Goddard Space Flight 
Genter or the EROS Data Center. Although the Landsat satellites theoretically 
make coverage of a particular area available every 9 days, the combination of the 
preceding two factors can easily limit availability to zero scenes. In the case 
of the Illinois lakes selected, although the Illinois EPA was interested in cover- 
age from June until September 1978, due to weather conditions, the periodic 
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shutdown of the landsat sensors and presumably data loss dui in# retrieval and 
recording* only two scenes were available over the selected lakes during this 
period. Of these scenes, neither was acceptable for processing as both had thin 
streams of small cumulus clouds running diagonally across the Images rendering 
them essentially useless. Owing to these problems, cooperation with the Illinois 
State EPA wa. necessarily lialted. 

Scene availability for the other areas of Interest was somewhat better. How- 
ever, coordination of ground sampling teams on selected flyover dates was plagued 
by weather, finances and lack of personnel available to sample. Table 1 1-9 pre- 
sents a list of Landsat scenes which were acquired for the study areas and Indi- 
cates the dates when ground truth sampling was performed on these lakes. 


2. The Ottawa-Nicolet Forest Lakes 

The U.S. Forest Service, Ottawa National Forest, had been monitoring the 
trophic status of lakes within the boundaries of their Jurisdiction. Interest in 
a landsat overview of their lakes was expressed in correspondence with USEPA per- 
sonnel at the Las Vegas Monitoring laboratory. Two scenes were acquired which 


Table 11-9. Landsat scenes acquired and sampling dates 


Area 

Scene I . D. 

Oa t e , 1978 

% Cloud 
cover 

Lake(s) 

obscured? 

Date of 
ground 
truth 

Kootenai National 
Forest and Flathead 
Lake (upper section) 

30146-17515 

7/29 

10% 

no 

8/07, 8/08 
8/25 

Flathead Lake 
(lower section) 

30146-17521 

7/29 

10% 

no 

8/07 

8/08 

Flathead lake 
(lower section) 

21293-17320 

8/07 

10% 

no 

8/07 

Nicolet National 
Forest 

21292-15432 

8/06 

10% 

no 

8/04 

Ottawa National 
Forest 

21293-15491 

8/07 

30% 

yes 

8/05 

Flaming Gorge Reservoir 

21288-17044 

8/02 

10% 

no 

8/02 

Flaming Gorge Reservoir 
(upper section) 

21287-16590 

8/01 

10% 

Part tally 

8/02 

Flaming Gorge Reservoir 
(lower section) 

21287-16592 

8/01 

20% 

Part ially 

8/02 
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covered the lakes of 111 r ' st In the Ottawa and Nlcolet National Foreat areas, 
loindsat 2 scene 21293-154 )i contained 30X cloud cover whlcii obscured most of the 
lakes selected for study. No further processing was performed on this area. 

landsat 2 scene 21292-134 32 contained no cloud cover over the area of In- 
terest. Processing was limited to color reconstruction of the Nlcolet National 
Forest urea. Figure 11-29 depicts the Nlcolet National Forest subsection from 
Land sat scene 21292-13432. Two o Inary masks were created using the program I.AKELOC 
to separate all the water bod i s In tills subsection from the surrounding land mass. 
Figure 11-30 and 11-31 are mask Images, In which the threshold for distinguishing 
water pixels was set at two different levels. As can be seen in these Images, the 
subjeit lakes were of such small size that further processing using only Landsat 
data would be of little use. Six bench mark lakes whlcii the Forest Service liad 
been monitoring were processed through 1.AKEL0C In preparation for classification 
and tiie results sent to USEPA in lais Vegas. It was determined tiiat due to the 
relatively small size of these lakes, unsupervised classification would not yield 
much useful information. No further processing was performed in this area. 


3. Flathead l.ake and the Kootenai National Forest 

USEPA personnel from the Western Energy office in Denver, Colorado, were 
planning to monitor Flathead lake In northern Montana during the summer of 1978. 

In addition the U.S. Forest Service, Kootenai National Forest, expressed interest 
in a Landsat overview of lakes located within their Jurisdiction. Accordingly 
efforts were made to coordinate sampling of these lakes with Landsat flyover 
dates. As can be seen in Table I 1-9, Flathead l.ake was sampled on August 7 and 8. 
The closest I*indsat scenes, covering the entire lake area, available for this time 
frame were acquired on July 29. As Flathead I.ake fell between landsat scanner 
positions it was necessary to acquire two consecutive scenes to reconstruct the 
entire lake surface area. A scene of the southern portion of Flathead Lake was 
also acquired which coincided with the August 7 sampling. However, the August 7 
landsat scene containing the northern portion of Flathead 1-ike was listP'* as un- 
available at the Goddard Space Flight Center. It is unknown whether thi was due 
to shutdown of the scanner or loss of data during processing by the center. 

Sixty-three lakes were Identified by the Forest Service as being of interest 
in the Kootenai National Forest region. Of these, nine were selected to serve as 
bench mark lakes for the study by agreement of EPALV and Forest Service personnel. 
The bench mark lakes were sampled by Forest Service personnel on August 25. Again 
the closest Iwindsat scene available covering these lakes was '.he scene acquired 
on July 29. All other scenes were listed as unavailable by the Goddard Space 
Flight Center. 

a. Pro cessing Kootenai National Fores t Lakes. Figure 11-32 shows a color 

reconstruction of the landsat 3 scene 30146-17515 for July 29, 1978. Flathead 
lake appears almost in its entirety in the lower-right-hand section of the image. 
The Kootenai National Forest area under study lies within the boundaries marked to 
the left of center within the image. Figure 11-33 shows the same area enlarged. 
Using LAKELOC, analysts at JPL located and isolated the 63 lakes of Interest. 

Figure 11-34 Is a binary mask of these lakes positioned exactly as they appear in 
the original Landsat scene. IPI. supplied EPALV with this linage and the statistical 
output from the prjgram STATUS for interpretation and comparison with ground 
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Figure 11-29. False color composite showing major 
portion of Nlcolet National Forest. 
Frame 21292-15432, August 8, 1978. 
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Figure 11-31. Binary mask of Htndsat MSS 

Channel 7. Second threshold. 
Land sat frame 21292-13432 
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Figure 11-32. False color reconstruction of Landsat frame 
30146-17515. Flathead lj»ke Is seen In the 
lower right corner. 


sampling measurements. Slr.ce ground measurements were acquired one month after 
the landsat flyover, analysis and comparison of these data were limited. Again 
due to the small size of the study lakes, no classifications were performed. 

b. The Fl athead L ake St udy Area. Flathead Lake Is the largest natural body 
of fresh water to be found west of the Mississippi River. The lake Is located in 
the northern portion of a large glacial valley near Glacier National Park In 
northwestern Montana. Flathead 1-ake Is 56.4 by 25.8 kilometers In '.ts extremes of 
length and width, and It has a maximum depth of i l 2 meters. (1) The lake shore- 
line Is 185 kilometers, more than half of which Is composed of rock and gn 
which extends lakeward to a depth of approximately 30 meters. (2) Appro*!- . iy 
50% of the shoreline Is composed of rocky cliffs. The steepness of much ( the 
drainage area and Its close proximity to the lake greatly limits Mie amount of 
drainage Into the lake. 

yi ORIGINAL PAGE 

ORIGINAL PAGE IS COLOR PHOTOGRAPH 

OF POOR QUALITY 



PRECEDING PAGE GLA 


MOT FILMED 


A portion of Landsat frame 30146-17513 showing 
the Kootenai Mlional Forest study area. 
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A binary mask of the 63 lakes of intere.it in 
the Kootenai National Forest 
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The Montana Water Pollution Control Authority has classified Flathead lake 
as an A-open Df lake (Ref. 1 1-15). This classification stipulates that the lake's 
water may be used tor drinking, culinary and food processing purposes suitable for 
use after simple disinfection treatment and removal of natural impurities. The 
water quality must also be maintained suitable for swimming and recreation. 

In terms of a temperature profile, Flathead lake has been classified as a 
cold dlmictlr lake which exhibits submerged depression individuality. This clas- 
sification Indicates that "there are various bays and depressions in the like each 
of which may liave Its own thermocllne that differs in position and thickness in 
the different depressions, and each depression may have Its own individual sea- 
sonal history. Furthermore, while much of the lake stratifies thermally, there 
are several extensive shallow bays which do not do so" (Ref. 11-15). 

c. Processing Flat he ad Lak e. As a response to a request by the EPA Western 
Energy Office for a demonstration of the capabilities of IPL's lake classification 
programs, a scene covering the northern half of Flathead lake was acquired for 
trial processing. The scene, 2879-17303, was recorded by landsat on June 19, 1977. 
This scene was logged into VICAR format using batch processing techniques which 
required approximately 25 minutes of computer run-time. On the following day the 
scene was processed through LAKELOC and classified with the programs USTATS and 
BAYES using interactive VICAR processing. The steps required 1.1 hours of computer 
run-time. The end product was a classification map reproduced in Fig. 11-35. 
Processing of this scene proved to the IPL's satisfaction that the lake classifi- 
cation programs could be used to produce digital classifications in a timely 
manner. However, this accomplishment acknowledged significant conditions which 
are listed as follows: 

(1) The subject lake should be of such dimensions as to warrant using 
digital classification techniques. 

(2) The time lapse between acquisition of data by Landsat and avail- 
ability to the user is such as to preclude any notion of "real-time" 
processing. 

(3) Such a classification is subject to interpretation based only on 
what can be discerned from the classification map in comparison 
with what is generally known about the lake's trophic status. 

There was no supportive ground truth available for Flathead Lake on June 19, 1977. 
Therefore this classification map represented an ex ellent overview of the lake, 
but the thematic classes represented are subject to interpretation without the 
benefit of supportive ground truth. 

Efforts were subsequently made to include Flathead Lake as a study area in 
the final demonstration phase of the Lake Classification projects. EPA personnel 
at the Western Energy Office in Denver agreed to provide bample data coinciding 
with landsat flyover dates. As indicated previously, Landsat data was limited to 
1 scene recorded on August 7, for which ground truth directly corresponded. 

Figure 11-36 shows a map of Flathead Lake with ground sampling sites as marked by 
EPA personnel. Two other scenes were also acquired which were recorded by Landsat 
9 days prior to the ground sampling. 
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Figure 11-35. An unsupervised classification map of Flathead Lake. Landsat frame 
30146-17515. 
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Figure 11-36. A conventional map of Flathead Lake. Water 


truth sites are marked on this map 
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Since Flathead lake ilea within two consecutive Land sat scenes. It was neces- 
sary to digitally connect two scenes to create one Image of the total lake surface 
area. Figure 11-37 Is a color composite Image constructed from the two Landsat 
scenes of Flathead lake acquired on July 29, 1978. The VICAR program STATS2 was 
used to gather mult lspe?tral statistics for each sampling site to be used In con- 
junction with analysis of classification statistics and ground measurements taken 
nine days prior to the landsat flyover. Table 11-10 lists the landsat statistics 
from STATS2. These products, accompanied by classification statistics, were sent 
to EPALV for analysis. As of this writing these data are still under interpreta- 
tion by EPA laboratory personnel. 

landsat scene 21293-17320 acquired on August 7, 1978, was also processed. 

As mentioned previously this scene contained only the lower half of the lake due 
to the landsat scanner position. Since corresponding ground truth measurements 
were available for this scene. It was classified using the unsupervised routine 
in USTATS. Three spectral classes were readily identified using USTATS. These 
mult Ispect ral statistics for each class were sent to EPALV for analysis. These 
data remain under study by EPA personnel. 


4. The Flaming Gorge Study Area 

The Flaming Gorge Reservoir is a long body of water located in Sweetwater 
Gounty, Wyoming, to the north and Doggett County, Utah, as its southern extreme. 
Flaming Gorge is approx imately 146 kilometers long with a surface area of 137.67 
km*-. The mean depth of the reservoir Is 33.9 meters with a maximum depth of 
153.0 meters (Ref. 11-16). 

EPA observation conducted at Flaming Gorge Reservoir indicate the following 
conclusions regarding its trophic condition. 

Flaming Gorge Reservoir ranked third in overall trophic quality 
among the 14 Wyoming lakes and reservoirs sampled in 197!) when com- 
pared using a combination of ?fx water quality parameters. 

5. Processing Flaming Gorge Reservoir 

Coordination of ground truth sampling and Landsat flyover was most success- 
ful for the Flaming Gorge Reservoir. Located on the border of Utah and Wyoming 
t lie reservoir is a significant site to warrant the use of Landsat data. Three 
Landsat scenes were acquired covering Flaming Gorge on August 1 and 2, 1978. The 
scenes were acquired or consecutive days due to the 14% overlap in Landsat 
scanner paths making available a "second look" on August 2. 

Two scenes were required to entirely cover the reservoir surface on August 1 
due once again to the Landsat scanner position. Figure 11-38 reproduces a com- 
posite image which was created from the two consecutive scenes acquired on August 
1. Cloud cover obscuc *d a portion of the reservoir surface. It was determined 
ttiat the presence of clouds and cloud shadow would affect the classification 
results. However, the decision was made to attempt an unsupervised classification 
despite the interference. Figure 11-39 presents a thematic map produced through 
USTATS and BAYES for the August 1 composite scene. As can be discerned, there 
was considerable confusion in class discrimination, which caused classes to 
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Figure 11-37. A color composite ImaRe of Flathead l^ke 
using two I^ndsat scenes (30146-1 7515 
and 30146-17521). 
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Table 11-10. 

landsat MSS 
sale (frames 

statistics from 
30146-17515 and 

training sites on 
30146-17521) . 

F la t head 

lake 

mo- 



Mean DN 




Sigma 



Training site 

Matrix size 

Green Red 

IR1 

1R2 

Green 

Red 

1 K 1 

IR2 

Ml.N 

5x5 

29.3 13.3 

5.2 

3.4 

2.1 

1.8 

1.1 

2.0 

ML 

5x5 

24.1 10.8 

5.7 

3.0 

1.6 

1.5 

1.3 

1.7 

YB 

5x5 

23.8 12.1 

6.8 

5.8 

1.4 

2.3 

0.9 

3.1 

BAB 

5 x 5 

25.8 14.4 

7.5 

6.1 

l.f 

1.5 

1.3 

2.0 

MLS 

5x5 

23.0 10.3 

5.0 

2.3 

1.3 

1.0 

1.8 

1.6 

PB 

5x5 

28.3 13.1 

6 . 6 

4.7 

2.6 

1.1 

1.2 

1.2 


become fragmented or speckled. This was probably due to atmospheric interference 
from shadows cast by clouds. The effects are most pronounced in areas of cloud 
concent rat ion . 

More reliable result s were obtained when scene 21288-17044, acquired on 
August 2, was processed. As shown in fig. 11-40, this scene was not affected by 
cloud cover. A» a result the mult ispectrai classification for this scene, re- 
produced in Fig. 11-41, does not exhibit the same confusion in classes. Class 
color distribution is fairly uniform throughout the reservoir. Corresponding 
ground truth measurements were available for this scene from EPA's Western Energy 
Office. In addition, the locations for sample sites indicated by EPA personnel 
were identified. Mult ispectrai statistics for each sample site were also 
acquired using STATS2. These stalstics are reproduced in Table 11-11. This in- 
formation was sent to EPALV for comment and analysis. As of this writing it 
remains utulur EPA study. 

The products presented in this section are examples of the results of digital 
image processing techniques which can be performed on water quality data supplied 
by Landsat. No attempt lias been made to analyze the accuracy of the multispectral 
classifications without expert input from EPA laboratory personnel. Until EPA 
analysis is complete no significant conclusions can be drawn regarding the pro- 
cessing performed during this study. 
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Figure 11-38. A color composite Image of Flaming Gorge 

Reservoir developed bv mosaicking two Landsat 
scenes together 1 2 1 .'*8 7 - 1 6590 and 2 128 7-1659?), 
August 1, 1978. 
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Figure 11-39. An unsupervised classification map of Flaming Gorge Reservoir for 
August 1, 1978. 
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Figure 11-40. A color composite Image of Flaming Gorge 


Reservoir, Landsat scene 21288-17044, 
August 2, 1978. 
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Figure 11-41. A supervised classifi- 
cation map o I Flaming 
Gorge Reservoir for 
August 2, 19'/ 8. 
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Table 11-11. Landxat MSS etatlatlca from training ultra on Flaming Gorge Krucr- 
volr (Frame 21288-17044) 


Training alte 

Matrix 

H lze 

Mean I)N 




Sigma 


Green 

Red 

I K 1 

IK2 

< .r 04-ii 

Red 

IR1 

IK.’ 

Site 

1 

2 

X 

2 

39.2 

37.2 

13.5 

3.0 

.8 

.8 

.8 

0 

Site 

2 

4 

X 

4 

31.0 

25.1 

10.1 

3.0 

l.n 

1.2 

1.8 

0 

Site 

3 

1 

X 

4 

43.0 

40.0 

15.0 

5.0 

1.8 

1.2 

1.2 

0 

Site 

4 

4 

X 

4 

23.3 

20.5 

8.4 

2.0 

1.1 

1.1 

1.6 

0 

Site 

5 

4 

X 

4 

2b. 2 

21.1 

8.8 

2.0 

1.0 

.8 

1.8 

0 

Site 

b 

5 

X 

3 

21.0 

15.6 

5.8 

1.0 

.9 

1.2 

1.2 

0 

Site 

7 

5 

X 

3 

20. b 

15.9 

6.5 

1.0 

1.3 

1.2 

.8 

0 

Site 

8 

3 

X 

3 

22.2 

lb. 4 

7.0 

2.0 

.b 

.8 

1.5 

0 

Site 

9 

5 

X 

3 

19.8 

13.3 

4.3 

1.0 

.9 

1.0 

1.0 

0 

Site 

10 

5 

X 

3 

22.3 

15.3 

5.8 

2.0 

.8 

.8 

1.2 

0 

Site 

11 

b 

X 

3 

21.0 

14. b 

16.2 

1.0 

.9 

.9 

1.7 

0 

Site 

12 

5 

X 

3 









Site 

13 

5 

X 

3 

21.7 

14.4 

6.9 

1.7 

1.4 

3.3 

5.7 

4.8 

Site 

14 

5 

X 

3 

22.5 

14.4 

4.8 

2.0 

1.3 

1.2 

1 .4 

1.0 

Site 

1 5 

3 

X 

3 

2) .9 

14.3 

7.1 

2.1 

1.2 

2.6 

5.7 

b.l 

Site 

1 

3 

X 

3 

21.7 

14.5 

3.3 

1.0 

.8 

1.2 

1.2 
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